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The importance of including colour octet contributions in describing decay and production of quarkonium states 
is briefly discussed for two cases : the radiative decays of the T and the production of J/ip in the inclusive B 
decays. It is shown how information on the non-perturbative matrix elements can be obtained by comparing the 
theoretical expressions computed at next-to-leading order in a a with the experimental data. 



1. Introduction 

It is now about four years that the non- 
relativistic QCD approach for the description of 
phenomena involving quarkonium production and 
decays, has been introduced PJ. Without any 
doubt, Bodwin, Braaten and Lepage have pro- 
vided us with a powerful and rigorous framework 
where calculations can be undertaken both at 
perturbative and non-perturbative level. Within 
this approach, the factorization between short- 
distance and long-distance physics allows to ex- 
press physical quantities as inclusive decay widths 
or cross sections as: 



Q(H) = ^c„(0„(i/)> 



(1) 



where in c„ are encoded the high-energy modes 
of the theory and so can be calculated perturba- 
tively in QCD, while the matrix elements (MEs) 
are non-perturbative and can be cither obtained 
by lattice simulation [|| or by comparison with 
experimental data. Moreover, by means of the 
power counting rules of NRQCD, it is possible to 
organize the MEs in a hierarchy based on the rel- 
ative velocity v <C 1 of the heavy-quarks in the 
bound state. In the end one is able to organize 
all these terms in a double expansion in q s and 
in v and to make predictions, in principle, at any 
given order of accuracy. 

With NRQCD in hands, much of the effort of 
the last years (for very nice reviews on more re- 
cent developments see for example ||]) has been 
devoted mainly to improve at NLO in a s the de- 
termination of the short distance coefficient c n for 
the most relevant terms in the v expansion and 



use the data available on quarkonium production 
and decays to extract information on the non- 
perturbative MEs. The aim of the game has been 
to check the cornerstone of the predictive power 
of the theory: the MEs should be insensitive to 
the details of the hard processes and depend in 
a universal way only on the bound state physics. 
This allows, in principle, to extract the values of 
the non-perturbative parameters from the exper- 
imental data on one process and then to use them 
to make predictions for others. 

In the following Sections I will discuss the ne- 
cessity of including the colour octet contributions 
at NLO to correctly predict and/or explain ex- 
perimental data for two cases of interest. 

2. The photon spectrum in bottomonium 
decay 

A consistent description of the photon energy 
spectrum in T — > 7 + X decay requires the in- 
clusion of the fragmentation components within 
the NRQCD factorization approach 0. The dif- 
ferential photon decay can be expressed in terms 
of a convolution between partonic kernels C a and 
the fragmentation functions D a ^ 7 : 



dr _ 

E 

o.=q,q-,g 



CJz) 



—C a (x, fi F )D a ^(-, fi F ) 



(2) 



1 Emission from final state light quarks has been surpris- 
ingly considered only recently by Catani and Hautmann 
in the CSM framework (and presented here in Montpellier 
in the QCD94 Conference ! [§) 
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Figure 1. Total colour octet contribution on 
the LO, colour singlet photon spectrum. Notice 
that neither NLO QCD nor relativistic effects are 
included in the singlet contribution. The overall 
normalization has been chosen as to reproduce 
the LO result of Ref. §. 



where z — E-y/rriQ is the rescaled energy of the 
photon ( rriQ is the heavy quark mass ) . The first 
term corresponds to what is usually called the 
'prompt' or 'direct' photon production where the 
photon is produced directly in the hard interac- 
tion while the second one corresponds to the long 
distance fragmentation process where one of the 
partons fragments and transfers a fraction of its 
momentum to the photon. The NRQCD expan- 
sion for the coefficients C{(x) reads: 



Ci 



Q 



(T|0(Q)|T) 



(3) 



where i = j,q,q,g and d[Q](x, a s (m Q ), ) are 
the perturbative coefficients. The NRQCD sum 
is performed over all the relevant Q states that 
contribute at a desired order v. The perturba- 
tive coefficients for the colour octet states are 



now known at NLO jjj. We can then investi- 
gate the phenomenological applications of colour 
octet states, if information on the NRQCD MEs 
can be obtained. To this aim we proceed in 
two steps: first we obtain estimates by solv- 
ing the renormalization group (RG) equations || 
and then we test the reliability of these esti- 
mates analyzing their impact on the observable 
i? M (T) = T(T -> had)/r(T -> m + aO, whose 
value is well-known both experimentally || and 
theoretically at NLO (both for singlet and octet 
contributions (lO) ) . As a result the RG estimates 
turn out to be smaller that one could expect form 
the naive velocity scaling rules and so legitimate 
doubts on their reliability arise. Nevertheless the 
above mentioned analysis shows that on one hand 
RG estimates are to be thought as lower limit 
while on the other, consistency between theory 
and experiment in total decay rates, strongly dis- 
favor much larger colour octet MEs. I remark 
here that these conclusions strongly rely on our 
present knowledge of the theoretical expression 
which is only at NLO. As a recent calculation 
shows lO], NNLO corrections can be very large. 
In fig. |l| are shown the various contributions to 
the spectrum. At low values of z the fragmenta- 
tion from octets is of the same order of magnitude 
of the LO colour singlet one. Contrary to LO ex- 
pectations in the framework of CSM ||, we con- 
clude that the decay of T into a photon would 
not be useful for estimating photon fragmenta- 
tion functions. For values of z near the end-point, 
breaking of the fixed-order calculation is manifest 
and the resummation of both short-distance co- 
efficient in a s and the non-perturbative MEs in 
v, is called for. Finally we notice that the over- 
all effect of octet states is at its minimum in the 
central region of the spectrum (exactly where the 
singlet LO direct contribution dominates) and so 
it should be used to make comparison with exper- 
imental data. Finally it has to be stressed that 
relativistic and higher order strong corrections to 
the singlet should be included to have a consistent 
theoretical picture at NLOR. 



2 The calculation of 
337 are in progress [1 
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3. The inclusive decay B — > J/ip + X 

In the case of a inclusive B decay into a char- 
monium state the following factorization formula 
holds @: 

T(B -> H+X) = ^2C(b^ cc[n]+x) (0 H [n}) ,(4) 



which is valid up to power corrections of order 
^QCD/mb,c- (To this accuracy it is justified to 
treat the B meson as a free b quark.) The pa- 
rameters (0 H [n\), defined in [jjj, describe the 
hadronization of a couple of heavy quarks into a 
charmonium state while the coefficient functions 
C{b — > cc[n] + x) describe the production of a cc 
configuration n at short distances and can be ex- 
panded in the strong coupling a s (fi) at a scale /x 
of order 2m c . The terms of interest in the AB = 1 
effective weak HamiltonianQ 



H eff = 



Of 
V2 



E 

q—s,d 



v; b v cq 



Cm A 



C[ 8 ] Og 



(5) 



contain the 'current-current' operators 

Ox = [c7 M (l-75)c][67 /i (l-75)9] (6) 
8 = [cT A 7m(1-75)c][6TV(1-75)9]. (7) 

A next-to-leading order calculation has been re- 
cently completed for all S and P states ]l3| but 
I will focus here only on the decay into a J/ip. 
At leading order in the velocity expansion the 
spin-triplet S-w&ve charmonium states are pro- 
duced directly from a cc pair with the same quan- 
tum numbers, i.e. n = ^S^. At order u 4 rela- 
tive to this colour singlet contribution, a ip can 
materialize through the colour octet cc states 
n = 3 S*f \ l S^\ 3 P^\ where the subscript 'J' im- 
plies a sum over J = 0,1,2. Note that even if 
the colour octet contributions arc suppressed by 
v 4 they must be included because the weak effec- 
tive Hamiltonian favours the production of colour 
octet cc pairs. The approximate relation holds 



C 2 



(8) 



which shows how the NRQCD suppression is com- 
pensated by the structure of the weak Hamilto- 
nian. In order to compare the final results with 

3 For simplicity I do not discuss here penguin contribu- 
tions, although they have been included in the final results. 



Table 1 

Extracted values for colour-octet MEs (in 
10 _2 GeV 3 ) for J/ip production (numbers in 
parenthesis are input values). 
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(1.06) 


3.1 


1.5 



experimental data, we need some input informa- 
tion on MEs. We use (Of ( 3 S"i)) = 1.16 GeV 3 
and the determination of (Of ( 3 <Si)) = 1.06 • 
10 _2 GeV 3 from direct J/ip production at large 
transverse momentum in pp collisions fll4j . The 
other two colour octet MEs are not yet well de- 
termined separately. Defining 



Mtesi s \ 3 pf) = (oteso))- 



we can reproduce the CLEO data Br (B — > J/ip - 
X) = (0.80 ±0.08)%, with : 

:(8) 3r>(8) 



Mf s (% \ 3 Pf) = 1-5±S:? • 10- 2 GeV 3 



(10) 

where a conservative choice for the errors has 
been made. It is interesting to compare the 
central values in ( pTj| ) and the upper limits 
with other determinations of the parameter 
M^^S^^P^). This is summarized in Tab. 
The central value in the first line is about a fac- 
tor 3 larger than the central values obtained in[l0[ 
As emphasized in Ref. Jl4| the Tevatron extrac- 
tion is very sensitive to various effects that af- 
fect the transverse momentum distribution. In- 
deed, Refs. |t5| quote smaller values compatible 
with, or smaller than the central values above. 
The total production cross section in fixed tar- 
get collisions probes Alf^S^, 3 Pj 8 ^) (assuming 
the validity of NRQCD factorization, which may 
be controversial). Given that a different com- 
bination of MEs enters, the values obtained in 
Ref. Jl6| are certainly consistent with the above 
central value. Considering the uncertainties in- 
volved in charmonium production in hadron col- 
lisions, the above upper limit on M3.1 is the most 
stringent one existing at present. 
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4. Conclusions 

We have shown, by discussing both a decay and 
a production process, that the effects of including 
colour octet states contributions are relevant. In 
particular we have given examples on how infor- 
mation on non-perturbative MEs can be obtained 
by the comparison between theoretical quantities 
(where short distance coefficients have been cal- 
culated at NLO) and available experimental data. 
This procedure has allowed us (a) to make predic- 
tions for photon energy spectrum in the T decays 
gaining information from the total and leptonic 
decay rates; (b) to give an estimate for some of 
colour octet MEs relevant for the J/tp produc- 
tion in B decays and so to check the "univer- 
sality" by comparison with the values extracted 
from other processes. As a final comment I would 
like to note that in both the above analysis we 
have found values of the MEs smaller than one 
could expect from the naive power counting of 
NRQCD. Whether this is a general feature of 
NLO calculations (as a preliminary analysis for 
fixed-target hadroproduction and photoproduc- 
tion indeed shows) and how this can be reconciled 
with the scaling rules is under study. 
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